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Gamow's explanation of the exponential decay law uses complex "eigenvalues" and exponentially 
growing "eigenfunctions" . This raises the question, how Gamow's description fits into the quantum 
mechanical description of nature, which is based on real eigenvalues and square integrable wave 
functions. Observing that the time evolution of any wave function is given by its expansion in 
generalized eigenfunctions, we shall answer this question in the most straightforward manner, which 
at the same time is accessible to graduate students and specialists. Moreover the presentation can 
well be used in physics lectures to students. 



I. INTRODUCTION 

In 1928 George Gamow considered the exponential decay of unstable atomic nuclei.^ His theo- 
retical description was based on solutions of the stationary Schrodinger equation 

d 2 

Hip = k -tp, where H = -—+V(x) (1) 

dx z 

and V is a compactly supported potentialP Gamow's key idea was to describe decay with eigen- 
functions G{x) that, asymptotically, behave like purely outgoing waves in the sense that 

lim (G(x) - e ±lzx ) = 0. (2) 

This idea lead him to "eigenfunctions" with complex "eigenvalues" z 2 = E — iT (E, T > 0) rather 
than real ones, i.e. HG(x) = (E — iT)G(x). Assuming that the Gamow function G(x) evolves 
according to the time dependent Schrodinger equation 

d 

i — G(t, x) = HG(t, x) = (E- iT)G(t, x), (3) 

it decays exponentially in time G(t,x) = e~ iEt *G(x). 

However, Gamow's description does not immediately connect with Quantum Mechanics. While 
Eq. |l]) appears there, too, in Quantum Mechanics eigenvalues are real and wave functions are 
square integrable. Gamow functions, on the other hand, belong to complex eigenvalues and are 
not square integrable. In fact, it is readily seen from their purely outgoing behavior (Eq. ([2])) and 
z = V E — iT having negative imaginary part, that Gamow functions have exponentially growing 
tails. Such a function is not square integrable. So the question is: How does Gamow's description 
of exponential decay connect with Quantum Mechanics? 

There are numerous mathematical articles concerned with this question, e.g. Refs. 3 6, From the 
articles it is, unfortunately, often not so easy to extract the clear and straightforward answer to that 
question. It is this: A Gamow function G{x) is approximately a quantum mechanical generalized 
eigenfunction (i.e. scattering state). Since generalized eigenfunctions govern the time evolution of 
square integrable wave functions which are orthogonal to all bound states, there are special initial 
wave functions, namely those which are approximated by a Gamow function and which therefore 
approximately undergo exponential decay in time. 
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Of course, this answer needs a bit of elaboration. We need to qualify the various "approxi- 
mations": First, generalized eigenfunctions do not have exponentially growing tails. Therefore, 
Gamow functions approximate generalized eigenfunctions only locally, e.g. on the support of the 
potential. The physical wave function, which undergoes approximate exponential decay must be 
square integrable and therefore can only be locally given by the Gamow function, too. Finally, 
approximate exponential decay in time means that neither for very small nor for very large times 
exponential decay holds. It only holds on an intermediate time regimeP 

Except for Ref. [HI the pedagogical accounts on Gamow's description of exponential decay we 
are aware of, usually only sketch its connection to the quantum mechanical description based on 
square integrable wave functions^]. The purpose of our note is to explain this connection in more 
detail. Compared to Ref. [HI which is a fairly complete discussion for a particular potential, we 
will stress the general principles underlying the connection between Gamow's description and the 
quantum mechanical description in a way which seems the most straightforward one, and which 
will be accessible to interested graduate students as well as specialists. The method presented here 
can well be taken as starting point for further results concerning for example higher dimensions. 
The presentation is also useful for teaching decay phenomena in physics courses. 

In the next Sec.[n]we give a heuristic explanation for Gamow decay, which will subsequently be 
rigorously discussed. At the end of Sec. [IT] we will give an illuminating example for our arguments. 



II. GAMOW FUNCTIONS AND THE TIME EVOLUTION OF SQUARE 
INTEGRABLE WAVE FUNCTIONS 



Gamow had the right intuition, "eigenfunctions" corresponding to complex "eigenvalues" do 
give rise to long lived square integrable states, which decay exponentially in time. However, their 
presence becomes only apparent in special physical situations. The prototype of a potential that 
creates such a situation is the double well potential (see Fig. [TJ . Wave functions initially localized 
inside the double well are long lived if the wells are high, because at potential steps they are partially 
transmitted and partially reflected; if the steps are high, reflection outweighs transmission. At each 
time of transmission, it is natural to view the transmitted portion as being proportional to what is 
left inside the double well and thus exponential decay appears naturally. 



X 
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Figure 1. Plot of the double well potential. 



This point of view shows that exponential decay is not at all a "tunneling" phenomenon, as 
it is often intuitively assumed. For a metastable state to occur, it suffices that a potential has 
steps at which a wave is partially reflected. The rectangular potential well Al^ is thus another 
example, which allows for unstable but long lived states that decay exponentially in time. Here 
le(x) denotes the indicator function on [—£,£}; le{x) equals one on [—£,£] and zero otherwise. If 
A is large, states initially localiz ed on top of the rectangular potential well will be metastable and 
will decay exponentially in timeP^ Gamow decay applies here as well and in Ref. [14] the decay 
is analyzed along the lines presented here. The picture of a wave being partially transmitted and 
partially reflected is somewhat hidden in Gamow's ansatz, but becomes more apparent when we 
relate the true quantum mechanical time evolution of the meta stable state to the Gamow function. 

If the double well (Fig. [TJ) allows for a Gamow function G, then the truncated version of it, 
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namely lgG yields a long lived square integrable initial wave function. The question we address is: 
Does lgG decay exponentially in time? In this section, we will explain that, on intermediate time 
scales, it actually does. At first we will consider general one dimensional potentials, but at the end 
(Section II C ) we will return to the concrete example of the double well potential. Throughout this 
note we will only consider potentials V, which have compact support contained in [— L, L] (note 
that L ^ £) and since exponential decay is a genuine scattering phenomenon, we also assume that 
the potentials have no bound states. 



A. Heuristic Argument 

Assume that V allows for a Gamow function G. We will give a two step argument, which shows 
that IlG decays exponentially in time when evolved according to the time dependent Schrodinger 
equation idtip — Hip. The general solution of the time dependent Schrodinger equation is given by 
e~ lHt tp. In the first step we will establish a generic connection between the time evolution of any 
square integrable wave function and the Gamow function G. In the second step we will use this 
connection to show that e~' lHt lLG decays exponentially in time. 

So, how does e~ lHt lLG evolve in time? To find an answer, we need a method that makes the 
time evolution palpable. For this purpose, we will use the method of expansions in generalized 
eigenfunctions <p(k, x), which applied to an arbitrary square integrable wave function ip yields 

/>oo 

tp(x) = / ip{k)<p(k,x) dk with (4) 



ip(k) = / ip{x)tp(k, x) dx. (5) 

J — oo 

These generalized eigenfunctions are bounded, but not square integrable solutions to the stationary 
Schrodinger equation 

H<p(k, x) = ( - + V(x)^j <p(k, x) = k 2 <p(k, x). (6) 
An expansion in terms of ip diagonalizes H in a completely analogous way as the Fourier transform 

.2 

diagonalizes — . The time evolved ip can thereby be expressed in a very concrete analytical way 



e 



im ip(x) = ^(k)(p(k,x)e- lk l dk. (7) 



Why should the time evolution expressed in terms of an expansion in generalized eigenfunc- 
tions ([7]) be related in any way to the Gamow function Gl Because both, the generalized eigen- 
functions ip as well as the Gamow function G, solve the stationary Schrodinger equation ([6|; the 
Gamow function for complex k 2 = E — iT with E, T > and the generalized eigenfunctions for real 
k 2 > 0. This suggests that in some sense ip ~ G, when the complex "eigenvalue" is close to the real 
axis (r 1). According to Eq. ^ generalized eigenfunctions behave like plane waves in regions 
where the potential is zero. Therefore, combining plane wave behavior and "near Gamow function 
behavior" , we make the ansatz 

^{k,x)^r](k)l L G(x) + e lkx . (8) 
We need to determine 77. Plugging Eq. Q into Eq. ([6]), we find 

■^+V\ Hk)l L G + e lkx ) w k 2 { V {k)l L G + e tkx ). (9) 
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Now, HIlG « z 2 1lG and --rie 1 ' 1 = k 2 e lkx , so we can rearrange the above equation, putting 



dx 

f{k,x) = V(x)e ikx , such that 



r,(k)l L G(x) « (10) 
Integrating both sides with respect to x, entails that 

where / is some analytic function and z 2 — E — iT. We find that the complex "eigenvalue" E — iT 
causes the generalized eigenfunctions cp(k,x) to have a pole, when continued to the complex k- 
plane. The Gamow function is the corresponding residue. This was the first step of our heuristic 
argument. 

In the second step, we will use Eq. (Ill in Eq. (l8|) to show that e~ %m \i J G decays exponentially 



in time. We only need to calculate the integral in the eigenfunction expansion Q. The heart of 
this calculation lies in the fact that the first summand on the right hand side of Q dominates when 
f< 1, because \rj(k)\ is much larger than \e lkx \ for k ps Rez. Therefore, 



<p(k,x)r*ri(k)l L G(x) (12) 



and hence 



l L G(fc)« / l L G(x)fj(k)l L G(x)dx^cfj{k), 

J —OO 

/CO 
f](k)r](k)l L G(x)e- lk2t dk 
-OO 



To solve the integral notice that it is essentially the Fourier transformation of the Breit-Wigner 
function l/[(fc 2 - E) 2 + V 2 ] . They differ only by the appearance of an additional function |/(fc)| 2 
and the fact that in Eq. ( 13 ) we integrate over k instead of k 2 . Therefore, we change the integration 



■5 



variable 



iHl 



l L G{x) « cl L G{x) J 
= cl L G(x) 
= cl L G(x) 



l/(*)l s 



\k 2 -(E-iT)\ 2 
„ \k 2 -(E-iT)\ 2 

\h+Vt)\ 2 „. 



dk 



e' lk 1 dk - 



de 



o (e-£) 2 +r 2 2^ Jo (e-E) 2 +T 2 ' 2^i 



|fc 2 -(£-ir)| 2 



e _i * * dk 



(14) 



We substituted fc = y/e in the first integral and k 



/e in the second. Due to the fact that the 



Breit-Wigner function is strongly peaked at e = E if T <C 1 (see Fig. Sj), the integrand in Eq. ( 14 ) 
is localized about E > 0. Hence, we can replace /(±i/e) and by their respective values at 

e = E, so that 



l L G(x) & cl L G{x) 



\f(+VE)\ 2 , l/(-V^)l 

2^ 



1 



(e - £) 2 + r 2 



rfc 



c'l L G(a;) e 



(15) 



where we have used that the Fourier transformation of the Breit-Wigner function is the exponential 
function. Thus, e~ lHt \[ J G decays exponentially in time whenever f< 1. 



B. Towards Rigor 



The crucial point in our heuristic argument is the form of the generalized eigenfunctions, i.e. the 
ansatz ^ together with the function r] from Eq. (11), because the time evolution is determined by 



the way the generalized eigenfunctions look like. How can this be made precise? 
First of all, the expansion in generalized eigenfunctions is defined by 



ip(x) = / ('i/j + (k)u + (k,x) + ip_(k)u-(k,x))dk with (16) 



4'±(k) — / ip(x)u±(k, x) dx, fce[0,oo) (17) 

J — oo 

and since it diagonalizes iJ, we have 

/>oo 

e- iHt i(j(x)= (^ + (k)u+(k,x)+^-(k)u-(k,x))e~ lk2t dk (18) 



for every square integrable wave function ip. Here u±(k,x) solve Eq. (|6|. Eq. (16 1 differs from 
Eq. Q by the appearance of an additional summand and by integrating from to oo rather than 
from — oo to +oo. This form of Eq. (16) is due to Schrodinger's equation ^ being an ordinary 
differential equation of second order. As such it has two linearly independent solutions for every 
k 2 G C. In particular, it has two linearly independent generalized eigenfunctions, u+{k,x) and 
ii-(k,x), for every k 2 > and both are needed to have a complete basisP^ For some potentials 
u+(— k,x) = u_(fc,a;), e.g. when V — for which u+ = e lkx /\/2iT and u_ = e~ tkx / ^/2ir . Then 
Eq. ( |16[ ) can be rewritten in shorter form (like Eq. Q) with an integral extending from — oo to +oo, 
but in general this is not possible. For more details on expansions in generalized eigenfunctions see 
Ref . M 

A potential with finite range perturbs the free Hamiltonian merely on its range, so that the 
generalized eigenfunctions are essentially plane waves, distorted only on the range of the potential. 
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Hence, we make the usuaf^ ansatz to obtain the precise generalized eigenf unctions: 

b+{k)e~ ikx 



a+(k)e lkx 
c + (k)e lkx 



f (k x) = l c -^ e ' ikx 



, x < —L 
, x > L 

, x < —L 
, x > L 



(19) 
(20) 



with k G [0, oo). Physically, Eq. (19) is a wave incident from the left together with its transmitted 
and reflected part; analogously Eq. (20) is a wave incident from the right. Now, it would be natural 
to set a± 
and \c±\ 

However, we choose a different route by setting c± — 1. Then the functions f±,a± and b± are 
analytic in k, because /+ (/-) is analytic for x > L (x < —L) and therefore needs to be analytic 
in k for all x. Thus f± have no poles. The poles occur upon normalization. Let u± denote the 
normalized eigenfunctions, then they need to satisfy 



1, thereby normalizing the amplitude of the incident wave to one, because then \b±\ 
12 would be the reflection- and transmission coefficients, respectively. 



which they do if 



S{k-k') 



u±(k,x) 



u±{k, x)u±(k' , x) dx, 



l_ f±(k,x) 
2tt a(k) ' 



k e [0,oo). 



(21) 



(22) 



Here a := a + = a_. To see that a+ is identical to a_, note that the Wronskian W(/+,/_) = 
f + f_ - f' + f- is independent of zP Thus, we can calculate the Wronskian once using the f± for 
x > L and another time using the f± for x < —L; the results must be identical and this leads us to 
a + = a_ . Note furthermore that the transmission- and reflection coefficients are now given by 



T{k) 



1 



1 



\a±(kW HkW 



and 



R±(k) 



\b±(k)\ 2 

Kfc)l 2 ' 



(23) 



respectively. Due to the fact that we set c± = 1, the analytic structure of the generalized eigen- 
functions is now evident: u± have poles whenever a — 0, because setting c± = 1 ensures that f± 
and a have no poles themselvesP 

We would like to argue now that Gamow functions are the residua of generalized eigenfunctions. 
From Eq. p2| we see that u± only have poles if a — 0. If z denotes a root of a(k), then the residue 
of u± at z is (up to a constant) given by f±(z, x). To see that f±(z, x) satisfies the purely outgoing 
boundary condition that Gamow functions satisfy (Eq. (|2|), observe that a — does not only 
cause the incident waves to vanish, but also b± — 1. For, W(f+, /_) = —2ika(k). So whenever a 
vanishes, the Wronskian vanishes, too. This is only possible if / + is a multiple of /_, which by a 
comparison of Eq. (19) with Eq. pO] ) implies that b±(z) = 1. Thus, f±(z,x) satisfies the purely 
outgoing boundary condition ([2|. Observing that f±(z,x) solves Eq. Q, too, we conclude that it 
is a Gamow function, so that 



Res(u±(k,x)) = cf±(z,x) = cG{x). (24) 

k — z 

Moreover, since 

HG{x) = Hf±(z,x) = z 2 f ± (z,x) = z 2 G(x), (25) 

z 1 is the corresponding eigenvalue. Eigenvalues corresponding to Gamow functions are usually 
called resonance. We will use this term not only for z 2 , but also for z. 
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Our ansatz (|8| is now rigorously replaced by the Laurent series of u±(k,x) about the root z of 
a(k) 

u±(k,x) = T}(k)G(x) + a„ (x) + af(x)(k - z) + ... , (26) 

where r)(k) = c/(k — z). But notice the difference to our ansatz (J8|, where the Gamow function 
comes with the factor 1/,, which truncates the exponential tails of G. In the Laurent Expansion (26 1 



this factor does not appear. Since u± is bounded and G increases exponentially, the remainder must 
compensate the exponential tails of G for large x, so that the whole sum remains bounded. The 
principal part of the Laurent Expansion therefore dominates only locally on bounded intervals. That 
is sufficient for our purposes, because we are interested in the exponential decay regime, where the 
majority of the wave function's mass lingers in intervals around the range of the potential. 

We wish to explain now under which conditions the principal part of the Laurent expansion 
governs the time evolution, so that exponential decay becomes dominant. For that purpose we 
consider the time evolution of ip = l^G, where G shall denote the Gamow function with eigenvalue 
z 1 = E — iT '. Separating the principal part from the remainder of the Laurent expansion (26 1, we 
obtain for the eigenfunction expansion ( 18 ) 



l L e- tMt l L G(x) =1 L 



(l L G + (k) + l L G_{k))n{k)G{x)e 

oo 

(l L G + (k)(u+(k,x) - v (k)G(x)) 



-ik A t 



dk 



(27) 
(28) 



+ l L G_(k)(u-(k,x) - r)(k)G(x)))e 



-ik t 



dk. 



Here we have inserted an additional factor 1^ for the reasons we have just explained. We have 
chosen the indicator function on the range of the potential only for simplicity. Regions growing 
with time would have been possible, tooPS 




Figure 3. The disk on which the Laurent series for u± converges. Note that the figure also illustrates that 
higher order resonances have larger decay rates. 



The Gamow contribution (27) dominates the time evolution whenever is strongly 

peaked at k = Re z and if T <C 1, because then (28) is negligible. We need both conditions, since 



the Laurent expansion of u± is only valid for those k for which the series converges. As V might 
have more than one resonance, the radius of convergence is in general finite, so that the Laurent 
expansion can only be used on a bounded interval 1 C 1 (see Fig. [3]). The contribution from 
this interval 1 to (28) can be neglected for small decay rates f « 1, since then the resonance is 



close to the real axis, in which case the remainder of the Laurent expansion is negligible compared 
to its principal part. The contribution from l c = [0, oo) \ I can be neglected when 1lG± is 
strongly peaked at k = Rez, because then 1lG ± is small on X c . That 1lG ± actually is peaked at 

k = Rez comes from the fact that 1lG± accumulates around k = Rez, because for these k we have 
u± w G when f<l. Putting the contributions from I and I c together, we see that neglecting the 



integral (28) causes a small error whenever f<C 1. 
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We refer to our heuristic Section II A to see now how the exponential decay arises from the 



principal part of the Laurent series. There is not much more to say on that. 



C. Example 

Now we consider the double well potential (see Fig. [I]). Also in this section we work with units 
in which to = 1/2 and h = 1 to keep formulas short. However, to provide a feeling for dimensions 
we will give some of the results in SI units. In this regard, recall that we study the dynamics 
of a-particles (He 4 , nuclei). The mass to that appears in Schrodinger's equation if it is written 
in SI units, therefore does not refer to electron mass, but the mass of a helium nucleus, which is 
6.69 x 10~ 27 kg. The conversion between our units and SI units then follows the scheme: xgi = a x, 
k S i = k/ao, A S i = |-Eu| A and t S i = (h/\Ei\)t, where a o = 7.2 x KT 15 m and E 1 = -0.1 MeV. 

Due to the simplicity of the double well potential, the generalized eigenfunctions can be calculated 
explicitly. For brevity we only give 

! [ (a(k) e lkx + b+{k) e~ ikx ,x<-(£+5) 
u + (k 7 x) — —j= — — - < ai(k) cos(kx) + i a^ik) sin(fca;) , |ac| < I (29) 



/27T a(k) 



e lkx , x > £ + 6, 



where a(k) = a\{k)a2{k) (the full formula and the auxiliary functions are given in the Appendix). 
This example perfectly illustrates the structure of the generalized eigenfunctions that we have found 



in Section II B In particular, it reflects that u + has poles whenever the function <x(fc) has roots. 
The decay rates (in SI units) of the corresponding metastable states are connected to the imaginary 
part of these roots by the formula 

rsi = T r = f Imz2 ' (30) 

where z denotes a root of a(k). 

Let us illustrate the decay rates of the double well potential using experimental data. Uranium 234 
for example, has an experimentally measured decay rate of r exp = 1.2936 x 10 -13 1/s (see Ref. \2U\ 
Tabic 1). To model a-decay of U 234 with the double well potential, we choose parameters which 
seem reasonable from what we know experimentally: For £ we choose the radius of the U 234 nucleus, 
i.e. £ = 1 (£ S i = 7.2 x KT 15 m). For S we choose the value 2£ (S S i = 14.4 x lCr 15 m). And A shall 
have the same order of magnitude as the Coulomb repulsion experienced by the a-particle at radius 
(■ (Vbouiomb(^) = 36 MeV), but its particular value is fitted such that the theoretical decay rate is 
close to the experimental decay rate. The best value is A = 436 (Agi = 43.6 MeV). We calculate the 
roots of a(k) and thereby the decay rates of the double well potential numerically (we have used 
Mathematica®). The root that describes the decay of U 23 best is 

hz sl = 1.0967 x 10~ 19 - % 8.5951 x 10 -55 (m/s) (31) 

and the corresponding decay rate calculates to Tgi = 1.3361 x 10~ 13 ^, which is very close to the 
experimental value. 



From the generalized eigenfunctions, we obtain the Gamow functions by Eq. (24 1 and the fact 



that b± = 1 if a = (see Sec. II B), 



~ izx ,x<-(£ + 6) 

G(x) — c ^ a±(z) cos(zx) + i 02(2) sin(zx) , |x| < i (32) 
y zx ,x>£ + 5, 

where z denotes a root of a(z). It is immediately evident that G satisfies the purely outgoing 
boundary condition pi). Moreover, we see that the Gamow functions have exponentially increasing 
tails, confirming our considerations in the introduction (Sec. |TJ). 
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Figure 4. The modulus of the generalized Fourier transform of leG (solid thin line) in comparison to |»j| 
(dashed thick line) is shown. For the figure we have used the same parameter values as for the Uranium 234 
example, which is discussed in the text (i.e. £si = 7.2 x 10~ 15 m, Ssi = 14.4 x 10 _ls m, Asi = 43.6 MeV). 



Here z denotes the root of a(k) from Eq. (31 1. We omitted the index SI for brevity. 



rj, which we have used in Sec. 
vanishes when evaluated at a resonance. 
liG(x) = clicos(zx) and 



II A 



via Eq. (17) 



Using Eq. (29) for the generalized eigenfunction and Eq. (32) for the Gamow function, we can 
calculate the generalized Fourier transfor m (lg G) 



This allows us to check that 
Due to the fact that a = aia%, either a\ or a 2 
Suppose z is a resonance for which a 2 vanishes, then 



(liG) + (k) = c I cos(zx) 



cos(fca;) 
a 2 (k) 



dx 



, . sin(fcx) 
cos(za;) _ —-. dx 



ai(fc) 



. .cosffca;) , 
cos(zx) — — dx 
a 2 (k) 



c ( sin((fc 
d 2 {k) \ k - 



z)t) 



sin((fc 
fc4 



z)t) 



(33) 
(34) 
(35) 



where we have used that cos(zx) sin(fcx) is antisymmetric in the second step. In Fig. [4] we have 
plotted the modulus of (1^G) + in comparison to the modulus of rj = c/(k — z). The amazing 
agreement confirms that (1^G)+ ~ rj. 



III. CONCLUSION 



Although, Gamow's approach based on complex "eigenvalues" directly leads to the exponential 
decay law, it involves exponentially growing "eigenfunctions" . Therefore, Gamow's approach to 
exponential decay can not be taken at face value. However, Gamow's exponentially growing "eigen- 
functions" are approximately quantum mechanical generalized eigenfunctions of the Hamiltonian. 
Since generalized eigenfunctions govern the time evolution of square integrable wave functions, 
Gamow functions give rise to special (square integrable) initial wave functions, which approxi- 
mately undergo exponential decay in time. So, Gamow's approach taken with a grain of salt yields 
in fact the explanation of the exponential decay law. 
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APPENDIX: GENERALIZED EIGENFUNCTIONS OF DOUBLE WELL POTENTIAL 



We now give the full formula for the generalized eigenfunctions of the double well potential 

'a(k)e ikx + b + (k)e~ lkx 
ci(fc) e lkx +c 2 (k) e~ lkx 
a\{k) cos(kx) + i a 2 (k) sm(kx) 
c 3 (k) e fkx + Ci (k) e~ fkx 



u+(k, x) 



1 1 

2tt o(fc) 



where 



a{k) = a\{k)a2{k) , A(k) = cos kS — - { — + — k5. 



x < -{1 + 5) 
-(£ + 5)<x< 
\x\ <£ 

£<x<£+5 
x> 1 + 5, 



2\k k 



a 1/2 (k) = e lM [A(k) + -e 



k k 



i fk 


k 


2\k 


k 



-ik2( 



ik2t\ 



kj\k k 



c 3/4 (fc) = Ui+^e^m+s) and ~ k =^Z7x. 



(36) 



(37) 
(38) 
(39) 

sin k5, (40) 
(41) 
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